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The critical nuclear charge Z c required for a heliumlike atom to have at least one bound state 
was recently determined with high accuracy from variational calculations. Analysis of the wave 
functions further suggested that the bound state changes smoothly into a shape resonance as Z 
crosses the critical value. Using variational calculations combined with the complex coordinate 
rotation method, we study the energy and width of the resonance for Z < Z 0 , thus providing direct 
evidence of the validity of this hypothesis. The variation of the resonance width with Z is found to 
be in good agreement with a model derived from analysis of the 1/Z perturbation series. 
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I. INTRODUCTION 


Let us consider the general problem of a two-electron atom with a nuclear 
described, in atomic units, by the Hamiltonian 


charge Z and infinite nuclear mass, 


H 0 = ~ 



2 ri r 2 r \2 ’ 


(1) 


where r* is the position of electron i with respect to the nucleus and r\ 2 the inter-electron distance. The scaling 
transformation —> Vi/Z , E —> E/Z 2 yields 


H = — 
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It is then possible to treat the electron-electron interaction as a perturbation, with a perturbation parameter A = 1/Z, 
leading to the following expansion for the ground-state energy E(Z): 


E(Z) = £ 


E n 


n =0 


(3) 


This perturbation series attracted considerable interest, being one of the few convergent perturbation expansions in 
atomic and molecular physics. A long-standing goal has been to determine the radius of convergence A* = 1/Z* as 
precisely as possible, and understand the relationship between Z* and the critical charge Z c , defined as the charge for 
which the ground state energy reaches the first detachment threshold, i.e. E(Z C ) = —1/2. An important step forward 
was the calculation of the first 402 coefficients of the perturbation expansion by Baker et al. [1J, whose subsequent 
analysis yielded the estimates Z* ~ 0.91103 [l| and Z* ~ 0.911028 26 [1|. On the other hand, a recent study of 
E(Z) using high-precision variational calculations 0 yielded Z c = 0.911028 224 077 255 73(4), later confirmed in [U, 
implying that Z* = Z c within the current accuracy of the determination of Z*. 

However, the nature of the transition occurringat the critical point has not been fully elucidated yet. Early studies 
by Stillinger 00 as well as more recent ones [7j connected with doubts on the accuracy of the determination of 
Z* [I] suggested that the system could persist as a bound state embedded in the continuum even for Z < Z c . This 
possibility was questioned by Reinhardt in an analysis based on the theory of dilation analyticity 0. In [3j an analysis 
of the variational wave functions showed that the outer electron remains localized near the nucleus even for Z cs Z c , 
and a simple model of the potential it experiences pointed to the existence of shape resonances. Overall the results 
of |3j are consistent with a smooth transition from a bound state for Z > Z c to a shape resonance for Z < Z c . 
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In this context, it is highly desirable to investigate numerically the existence and properties of resonances for 
Z < Z c . Resonances were studied with the complex rotation method in 53. but in that work only the region 
A > 1.11 (i.e. Z < 0.9009) was explored, still relatively far from the critical value. Here we present more accurate and 
complete results where the range of investigation was extended up to Z = 0.9103. Our results are in full agreement 
the hypothesis of a smooth transition from a bound state to a resonance formulated in I, and put much stronger 
constraints on the possible existence of a bound state for Z < Z c . Furthermore, a model of the resonance width 
derived solely from the analysis of the asymptotic form of the 1/Z perturbation coefficients na is in good agreement 
with our data, and gives further indication that the state of the system starts acquiring a finite width exactly at 

^ = Ze¬ 


ll. NUMERICAL CALCULATION OF RESONANCES 


In order to obtain the energies and widths of resonances, we diagonalize numerically the complex rotated Hamilto¬ 
nian m in a variational basis set. We use the perimetric coordinates (denoted x , y, z) and Sturmian wave functions 


A- 


a,(3 


^x,y,z) = (-1) 


— \ n x+n y +n z 
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(4) 


where n x ,n y ,n z are non-negative integers and L n the Laguerre polynomials fl2j . The variational basis is truncated 
by imposing the conditions 


n x + n v + n z < n r> 


n x < n r ‘ 


(5) 


This approach allows to write the Hamiltonian in the form of a sparse-band matrix, and gives high numerical 
accuracy with double-precision arithmetic 5153- It has been shown to yield very precise results for the complex 
energy of resonances [15j. Finally it is well suited for studies where parameters of the physical system are changed 
(such as, here, the nuclear charge Z). Indeed, the small number of variational parameters (the two length scales cU 1 
and /? _1 ) makes it easy to readjust them for each value of Z. 

Our results are summarized in Table U and Figure [1] We cross-checked the determination of Z c in Refs. II by 
computing the energy for Z = 0.911028 244 077 2 and Z = 0.911028 244 077 3 (see the last two lines of Table |I|. In 
both cases all 15 digits are in agreement with those of [!]. We thus confirm that Z c lies between these two values; 
a more detailed analysis is possible but would lie outside the scope of the present work. This also shows that our 
method can reach an accuracy of about 1 x 10 -15 a.u. that is essentially limited by numerical noise. 

In the region Z < Z c , we checked the closest value to Z c studied so far 0, i.e. Z = 1/1.11 = 0.9009009... and 
find 


Re(£’) = -0.497131212 367 59, Im (E) = -4.997 292 209 10“ 5 . (6) 

This is in agreement with, and much more accurate than the result of [lOj, Im(F') = —6(3) 10 -5 . The study was then 
pursued by increasing Z in small steps, going as close to Z = Z c as possible. As can be seen from Tableland FigureQ] 
the width of the resonance drops rapidly with increasing Z. The ~ 10~ 15 a.u. numerical accuracy of the calculation 
implies that only resonance widths in the 10~ 14 a.u. range or larger can be obtained with sufficient accuracy. This is 
what limits the explored range to Z = 0.9103 where the width is 4.7(1) 10 -14 . At Z = 0.9104 the width is already 
too small to be resolved. 

Overall, our results confirm that the bound states is transformed into a resonance. In addition, the parameters a, /3 
of the variational wave function are found to vary smoothly as a function of Z 1 in agreement with what was observed 
in I, suggesting a smooth transition from bound to resonant state. Rigorously speaking, there is no proof that this 
transition occurs exactly at Z = Z c , since the interval Z £ [0.9103, Z c ) was not studied. One could improve further 
the accuracy of the calculation by using quadruple precision arithmetic and a larger variational basis. However, in 
view of the extremely fast decrease of Im(A), this would not allow to get much closer to Z c and would not be worth 
the increased computational effort. Whatever the accuracy of the calculation, some small range in the vicinity of the 
critical point would remain beyond reach. Thus, the only clear-cut conclusion is that 

Z* > 0.9103 (7) 

However, there is no physical argument to support the hypothesis that, for decreasing Z. a bound state would 
persist for some small range [Z *, Z c ] and then disappear at Z = Z *. As discussed in jS, [l6[, such a scenario is 
compatible with the theory of dilation analyticity only if this bound state is not square-integrable. This possibility 
seems to be ruled out by the observed smooth variation of the parameters describing the wave function. 

In order to gain further insight on the nature of the atomic state in the range Z £ [0.9103, Z c ] and estimate the 
value of Z at which the width of the resonance goes to zero, we will use a model of the resonance width derived by 
Ivanov and Dubau 52- This is the object of the next Section. 
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z 

R e(E) Im (E) 


0.905 

-0.498 261 103 900 07 -6.224 027 27 

10 “ ub 

0.906 

-0.498 542 502162 97 -2.552 367 30 

l 0 -°6 

0.907 

-0.498 826 879 424 93 -7.482 477 6 

IQ-07 

0.908 

-0.499 114 312122 79 -1.203 842 0 

IQ-07 

0.908 5 

-0.499 259111008 48 -3.233 035 

10-08 

0.909 

-0.499 404 559 766 78 -5.397 92 

Iq-09 

0.909 5 

-0.499 550 593 342 34 -3.893 5 

10 -io 

0.909 8 

-0.499 638 468 696 42 -3.850 

10- 11 

0.910 

-0.499 69715161159 -4.79 

10~ 12 

0.9101 

-0.499 726 521799 42 -1.32 

10 -12 

0.910 2 

-0.499 755 910 711 12 -2.9 

10- 13 

0.910 3 

-0.499 785 318 046 200(1) -4.7 (1) 

10- 14 

0.911028 224 077 2 

-0.499 999 999 999 983(1) - 


0.911028 224 077 3 

-0.500 000 000 000 013(1) - 



TABLE I: Real and imaginary part of the energy of the lowest resonance of a two-electron atom of nuclear charge Z (in atomic 
units). Typical values of the parameters are: n rnax = 160, n™ ax = 25 (yielding a basis size N = 145 639), a ~ 1.5 — 3, 
P ~ 0.4 — 0.8, complex rotation angle 9 ~ 0.3 — 0.6. All digits are significant unless otherwise noted. A larger basis set (up to 
n max _ 2 qq^ n ™ ax — 33 ^ yielding N = 293 301) was used for the last three lines in order to get one more significant digit. 


Im(E) 



FIG. 1: Imaginary part of the energy of the lowest resonance of a two-electron atom as a function of the nuclear charge Z. The 
solid line is a fit by expression (Hill (see text for details). 


III. MODEL OF RESONANCE WIDTH 

Let us briefly summarize the approach of Ivanov and Dubau [ill ] . They first derived a dispersion relation connecting 
the real and imaginary parts of the E(Z) function, relying on the following analytic properties: 

(i) E(Z) is a regular function for Z > Z* (since the series ([3]) is convergent). 

(ii) It has an essential singularity at the point Z*. An indication of this behavior was obtained in [l] from approximating 
E( A) by a known function whose Taylor coefficients match very closely the perturbation coefficients. 

(iii) It has a second-order pole at Zq = 0, as shown in E3- „ 

(iv) It has a third singular point on the real axis for some value Z 2 , which was estimated to Z^ ~ 0.11 in El- 
The dispersion relation reads 


C-! , C - 2 


E{Z) = E 0 + —^ + 


Z 2 


— lim 
7r £->0+ 


r Z * 


lmE(t — ie) 


dt, 


z 2 


Z-t 


( 8 ) 
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where C-% and C _2 are real coefficients. From this expression one can deduce the following formula for the coefficients 
of the expansion 0: 


1 f z * 

E n = — lim Im E(t — ie ) f™ -1 dt 
tt Jz 2 e->0+ 


The following asymptotic large-n behavior of the coefficients E n was established in [J]: 


(9) 


E n ~ C {Z*) n n b e~ a ^ (1 + (10) 

V n / 

Ivanov and Dubau showed that the leading term of this asymptotic law can be reproduced if one assumes that the 
behavior of the imaginary part of the energy near the point Z* (for Z < Z*) is given by: 


Im(£ l ) ~ A 



( 11 ) 


where 


3 a 2 

p = —2 b — — and c = —. (12) 

We performed a least-squares fit of the coefficients E n from [I] in the range n G [101,401] by the first term in 
expression m (more precisely, we fitted ln(£'„) by the logarithm of that term with equal weights for all data points); 
Z* was set equal to 0.911 028 22. We find 


a = 0.26374(5), b= -1.9896(4). (13) 

These values are close to those obtained in 1| using a different procedure, 

a « 0.272 , b « —1.94. (14) 

Then we performed a least-squares fit oi\m[E{Z)\ by (fill) in the range Z G [0.905,0.9103] (again, we fit ln(Im[S(Z)]) 
by the logarithm of (1111) 1. Here, Z* was taken as a fitting parameter. This yields 

p = 2.42(9), c = 0.0180(5) and Z* = 0.911 276(12). (15) 

The values of p and c are in very good agreement with those deduced from (fl2l) and (fl3|l . p = 2.48 and c = 0.0174 
(or p = 2.38 and c = 0.0185 if one adopts the values of a and b from (flTl) reported in jlj). The value of Z* deduced 
from this fit is slightly too large, by 2.5 10 -4 (Z* cannot be larger than Z c ). We also tried fitting by a refined 
version of ED which was derived in Ref. 0 to reproduce also the next-to-leading-order term of ED. but this did 
not improve the agreement. This discrepancy could be related to the numerical difficulties discovered by Baker et 
al. [I]]: the asymptotic behavior of the perturbation coefficients E n manifests itself only at very high orders, and the 
nature of the singularity of E(Z) is revealed only extremely close to Z = Z*. It is thus possible that a more detailed 
investigation of the asymptotic behavior of E n could lead to a more accurate model of the variation of Im(£J) vs. Z. 

Despite the discrepancy on Z* , the good quality of the fit suggests that the expression ED accounts for the essential 
features of the variation of Im(£'). In particular, it shows that Im^E 1 ) does not drop to zero immediately beyond the 
last investigated value at Z = 0.9103 but continues to have a nonzero values up to a close vicinity of Z c . This 
strengthens the hypothesis of the ground state acquiring a finite width exactly at Z = Z c , implying the identity 
Z*=Z C . 


IV. CONCLUSION 

We have studied the atomic state in two-electron atoms slightly below the critical nuclear charge Z c = 0.911 028... 
by variational calculations combined with the complex rotation method. It was shown that already in a close vicinity 
of Z = Z c , the bound state is transformed into a resonance whose width increases very rapidly as Z is decreased 
further away from Z c . The variation of the width as a function of Z is in good agreement with a model derived 
from the analysis of the 1/Z perturbation series, suggesting that the state acquires a finite width exactly at Z = Z c . 
It would be interesting to get an independent confirmation by a recalculation and analysis of the 1 /Z perturbation 
series, up to even higher orders than what was done in [l|. 
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